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Matter wave switching in Bose-Einstein condensates via intensity 
redistribution soliton interactions 

S. Rajendran, 1 P. Muruganandam, 2 and M. Lakshmanan 1 

x ' Centre for Nonlinear Dynamics, Bharathidasan University, Tiruchirapalli 620024, 
India 

^Department of Physics, Bharathidasan University, Tiruchirapalli 620024, India 

Using time dependent nonlinear (s-wave scattering length) coupling between the components of a weakly 
interacting two component Bose-Einstein condensate (BEC), we show the possibility of matter wave switching 
(fraction of atoms transfer) between the components via shape changing/intensity redistribution (matter 
redistribution) soliton interactions. We investigate the exact bright-bright iV-soliton solution of an effective 
one-dimensional (ID) two component BEC by suitably tailoring the trap potential, atomic scattering length 
and atom gain or loss. In particular, we show that the effective ID coupled Gross-Pitaevskii (GP) equations 
with time dependent parameters can be transformed into the well known completely intcgrablc Manakov 
model described by coupled nonlinear Schrddingcr (CNLS) equations by effecting a change of variables of 
the coordinates and the wave functions under certain conditions related to the time dependent parameters. 
We obtain the one-soliton solution and demonstrate the shape changing/matter redistribution interactions 
of two and three soliton solutions for the time independent expulsive harmonic trap potential, periodically 
modulated harmonic trap potential and kink-like modulated harmonic trap potential. The standard elastic 
collision of solitons occur only for a specific choice of soliton parameters. 

PACS numbers: 03.75.Mn, 03.75.Lm, 05.45.-a 



I. INTRODUCTION 

The past decade has witnessed a considerable increase 
of interest in the experimental and theoretical studies of 
matter wave solitons of the darker— and bright^ 5 - types 
in Bose-Einstein Condensates (BECs). These solitons 
have attracted a great deal of attention in connection 
with the dynamics of nonlinear matter waves, includ- 
ing soliton propagation^, vortex dynamics^, interfer- 
ence patterns^ and domain walls in binary BEC o 10 ' 11 . 
From an experimental BEC point of view, bright solitons 
are created themselves as condensates^ 5 - while dark soli- 
tons exist as notches or holes within the condensates^. 
Note that bright solitons propagate over much larger dis- 
tances than dark solitons. The bright matter-wave soli- 
ton trains (multi-solitons) were experimentally observed 
by Strecker et al4 in 7 Li and Khaykovich et al£ in 87 Rb 
by magnetically tuning the atomic scattering length from 
repulsive to attractive nature, through Fcshbach reso- 
nance^. The recent experiments at Heidelberg and Ham- 
burg universities have shown the formation of dark soli- 
tons, their oscillations and interaction in single compo- 
nent BECs of 87 Rb atoms with confining harmonic po- 
tentiaU2rJ£. 

Many aspects of the above novel and experimentally 
accessible form of matter have been since then inten- 
sively studied; one of them concerns with the investiga- 
tion of the behaviour of multi-component BECs, which 
have been experimentally studied in either mixtures of 
different hypcrfine states of the same atomic species or 
even in mixtures of different atomic species. Experimen- 
tal generation of two-component BECs of different hy- 
pcrfine states of rubidium atoms in a magnetic trap^ 
and of sodium atoms in an optical trapi 7 - stimulated 



theoretical studies devoted to the mean-field dynamics 
of multi-component condensates^. Recently, Zhang et 
al.— proposed a method for independent tuning of scat- 
tering lengths in multi-component BECs. When a con- 
densate is cigar shaped and has relatively low density, 
that is, when the healing length of the components is 
much larger than the transverse dimension of the con- 
densate and much less than its longitudinal dimension, 
the transverse atomic distribution is well approximated 
by the Gaussian ground state and the system of coupled 
ID Gross-Pitaevskii (GP) equations is adequate to de- 
scribe multi-component condensates. 

In the recent literature, there has been a growing 
interest, both from experimental as well as theoreti- 
cal perspectives^— , in studying the dynamics of two- 
component BECs coupled to the environment such as 
external thermal clouds which leads to the mechanism of 
loading (gain) external atoms (thermal clouds) into the 
BECs by optical pumping or continuously depicting (loss 
of) atoms. It is interesting to note that matter wave soli- 
tons in multi-component BECs hold promise for a num- 
ber of applications, including the multi-channel signals 
and their switching, coherent storage and processing of 
optical fieldaSSr^S. 

Further, there has been increased interest in recent 
times in studying the properties of BECs with time vary- 
ing control parameters, such as (i) the temporal varia- 
tion of atomic scattering length which can be achieved 
through Fcshbach resonance 2 ^— , (ii) inclusion of appro- 
priate time dependent gain or loss terms^ -25 ' 33 " — , (iii) 
the temporal modulation of trap frequencie s) 25 ' 38 " — and 
so on. In particular, the study of matter wave solitons 
under time varying control parameters is one of the cur- 
rent active research field o 25 ' 38 ' 42 " — . Being motivated by 
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the above considerations, in this paper, we study the dy- 
namics of the exact bright-bright matter wave solitons 
of two component BECs under the above time varying 
control parameters. 

Most of the theoretical studies on the matter wave soli- 
tons in multi-component BECs have been carried out 
using cither numerical or approximation methods. For 
example, two component dark-bright solitons have been 
reported in^; bound dark solitons have been numeri- 
cally studied in^S, where it has been found that the cre- 
ation of slowly moving objects is possible; a diversity 
of other bound states has been generated numerically 
in^i. The present authors have obtained the exact dark- 
bright soliton solutions and their different kinds of in- 
teractions of the two component BECs by suitably tai- 
loring the trap potential, atomic scattering length and 
atom gain/loss^. We also point out that Babarro et al£2- 
have shown the possibility of switching phenomenon of 
matter- wave solitons via bright-bright soliton interaction 
of different species of two component BECs essentially 
without potential term or any modulation of control pa- 
rameters, through an approximate theoretical analysis 
and numerical simulations. In the present paper, we 
have investigated a different kind of matter wave switch- 
ing phenomenon (in two different hypcrfine states of the 
same species such as 85 Rb with equal intra- and inter- 
species atomic interactions) via intensity redistribution 
of exact bright-bright soliton interactions with the tem- 
porally modulated control parameters analytically. 

Specifically, in the present paper we bring out the exact 
bright-bright one-, two-, three- and ./V-soliton solutions in 
two-component BECs by simply mapping a class of two 
coupled effective ID CP equations onto the completely 
integrable two coupled nonlinear Schrddingcr (2CNLS) 
equations (Manakov system) and making use of the ex- 
act /V-soliton solutions of the latter system. In particu- 
lar we have demonstrated the bright-bright shape chang- 
ing/matter redistribution of two and three solitons un- 
der collision, while elastic collision occurs for a very spe- 
cial choice parameters. We have also shown the shape 
restoring property in the case of three soliton interac- 
tion. These types of elastic and shape changing inter- 
actions of two and three solitons have been well studied 
in the context of optical computing, where the intensity 
of light pulses are transformable between two modes of 
Manakov type systems^—. From the BEC point of view 
the shape changing interaction can be interpreted as the 
transformation of the fraction of atoms between the com- 
ponents, which is the so called matter wave switch. Such 
matter wave switching phenomenon can be used to ma- 
nipulate matter wave devices such as switches, logic gates 
and atom-chip2£r— . One of the long term prospectives 
of matter wave devices is their potential application to 
quantum information processing, for details see for exam- 
ple Refs.— ~— . In the present paper we have shown that 
such a matter wave switching phenomenon in two com- 
ponent BECs is possible via shape changing soliton in- 
teractions by suitably tailoring the trap potential, atomic 



scattering length and gain or loss term. 

This paper is organized as follows. In Section [TTl we 
present the ansatz for the two coupled GP equations in 
ID to be mapped onto the integrable 2CNLS equations. 
In Section IIIII we deduce the one-, two-, three- and N- 
soliton solutions of the two coupled GP equation from 
the one-, two-, three- and ./V-soliton solutions, respec- 
tively, of the 2CNLS equations. In section ITVl we bring 
out the one soliton solution, elastic and shape changing 
soliton interaction of two and three soliton solutions for 
different forms of the trap potential, gain/loss term and 
interatomic interaction (scattering length) . We have also 
shown the shape restoring property in three soliton inter- 
actions. Elastic collision occurs only for a specific choice 
of parameters. The analysis can also be extented to N- 
solitons without much difficulty Finally, in Section [Vj 
we present a brief summary of our study. 



II. ANSATZ FOR MAPPING TWO COUPLED GP 
EQUATIONS ONTO MANAKOV SYSTEM 

We consider the dynamics of a ID two-component 
trapped BEC with gain or loss term by the mean-field 
equations for the wave functions, "fi and ^2, of the con- 
densates II) and 12): 



1 d 2 V 3 

2 dx 2 



+ V(x,t)-^ + - 1 (t) 



J = 1,2, 



k=l 



(1) 



where V(x,t) = n 2 (t)x 2 /2 is the external time varying 
trap potential, which is expulsive for Q. 2 (t) < and con- 



fining for n 2 (t) > 0. Here n 2 (t) = u 2 x {t)/u] 



is the temporally modulated axial trap frequency, uj± 
is the time independent radial trap frequency, R(t) = 
2a s {t) / as, a- s (t) is the magnitude of the s-wave atomic 
scattering length, ag is the Bohr radius, Ujfc's are the 
signs of the s-wave scattering lengths, which are nega- 
tive for attractive and positive for repulsive interactions, 
fij is the chemical potential of the jth component and 
7(i) = T(t)/uj±, where T(t) is the gain (for +ve) or loss 
(for — vc) term, which is the phenomenologically incor- 
porated interaction of external atoms (thermal clouds). 
The time dependent gain/loss term corresponds to the 
mechanism of continuously loading external atoms into 
the BEC (gain) by optical pumping or continuous deple- 
tion (loss) of atoms from the BEC22~— . 

The atomic scattering length of alkali atoms such as 7 Li 
and 85 Rb atoms can be experimentally varied by suitably 
tuning the external magnetic field through the Fcshbach 
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resonance as 



a a (t) 



B(t) - B c 



(2) 
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where <x° is the scattering length of the condensed atoms, 
B(t) is the external time varying magnetic field, B is the 
resonance magnetic field and A is the resonance width. 
Recently, Zhang et al* 1 ^ proposed a method for a simi- 
lar kind tuning of scattering lengths for two component 
condensates of two different hyperfine states of 87 Rb. 
In the present study we have considered the attractive- 
attractive {(Jjk = — 1) two component BECs as in the 
case of two different hyperfine states of 85 Rb with equal 
intra- and inter-species atomic interactions. 
By applying the following transformation, 



^j(x, t) = exp 



i(ijt+ J ^Y~dt 



Qj(x,t), (3) 



where j = 1,2, Eq. (J]) can be transformed to a system 
of two coupled GP equations without gain/loss term and 
chemical potential as 



R(t)Y,^k\Qk\ 2 + V(x,t)- 



k=l 

where j = 1 , 2 and 



1 d 2 
2dx^ 



Qh (4) 



In order to solve for the unknown functions A, F and G 
in the above equations (|10[) we assume the polar form 



A = r(x,t)exp[i9(x,t)]. 



(11) 



One can immediately check from the relations (|10c[) that 
r is a function of t only, r = r(t), since G and R are 
functions of t only. Then from Eqs. (fTUf one can easily 
deduce the transformation function A given by and 
the transformations G(t) and F(z,t) through the follow- 
ing relations, 

r 2 = 2r 2 R, (12a) 

9 = -^tx 2 + 2br\Rx - 2b 2 r\ [ R 2 dt, (12b) 
2R J 

F(x, t) = X = V2r Rx - 2\/26tq / R 2 dt, (12c) 



G{t) =T = r\ 



R 2 dt. 



(12d) 



Here b, rg are arbitrary constants, and R and Q 2 should 
be related by the following condition, 



R(t) = exp 



j(t)dt 



R(t). 



(5) 



Eq. (|4]) can be mapped onto the 2CNLS (Manakov) equa- 
tions under the following transformatio n 38 i 43 i 48 i 59 i 60 , 



i (x,t) = Aq J (X,T), i = 1,2, 



(6) 



where the new independent variables T and X are chosen 
as functions of the old independent variables t and x as 



T=G(t), X = F(x,t), 



(7) 



while A = A(x, t) is a function of t and x. Applying the 
above transformation ([5]), so that 



d_ 
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dx 2 ~ xx dX 
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F 



x dx 2 



^dX 1 
_dG dF 



dt 



dx 



one can reduce Eq. Q to the 2CNLS equation of the 
form 




(13) 



which is a Riccati type equation for R t /R- Eq. ^ is 
the celebrated Manakov system^— , which is well stud- 
ied in the context of nonlinear optics, biophysics, plasma 
physics etc. Eq. (J9j is a completely integrable soliton 
system and exhibits interesting one-, two-, three- and 
A^-soliton solutions of bright-bright typ o 53 i 55 ~ — . From 
the solutions of Eq. ((9]), one can straightforwardly con- 
struct the one-, two-, three- and N- bright-bright soliton 
solutions for Eq. (TJJ), provided R t , 7(t) and satisfy 
Eq. ([T3]) . In the context of BECs, it is of fundamen- 
tal interest to study the bright-bright soliton solutions 
of Eq. (TTJ) . In the following we shall describe the multi- 
soliton solutions corresponding to the coupled GP equa- 
tion (Hj. 

Note that in Eq. (fTJ) , the variable x actually represents 
x/a±, where a± = h/(mu>±), and similarly t stands for 




qj =0, .7 = 1,2, (9) 



III. MULTI-SOLITON SOLUTIONS OF TWO COUPLED 
NONLINEAR SCHRODINGER EQUATIONS 



subject to the conditions that the functions A, F, G, il 
and R should satisfy the following set of equations, 



iA t + 
iAF t 



2 Axx 



^x 2 A = 0, 



(2A.F, + AF XX ) = 0, 



F 2 f? 
G t = =-2- = -|A| 2 . 
2 2 1 1 



(10a) 
(10b) 

(10c) 



The 2CNLS equation © in contrast to the sin- 
gle component NLS system admits solutions which ex- 
hibit certain novel energy sharing (shape changing) col- 
lision a 53 ' 55 " — . Recently, the general expression for N- 
soliton solution of the Manakov system in the Gram de- 
terminant form has been given in Ref^i by using Hirota's 
bilinear method. 

In order to write down the multi-soliton (Af-soliton) 
solution of the focusing Manakov system we define 
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the following (1 x TV) row matrix C s , s = 1,2, (2 x 1) 
column matrices ipj, and (TV x 1) column matrix 0, j = 
1,2,..., N, and the (TV x TV) identity matrix /: 



C, 



(4 



s) (s) 
,<*2 % 



= (0,0,...,0), 



(14a) 



with j = 1,2. The explicit form of the two-soliton solu- 
tion can be written as 



where 



i=i 







(20) 







(1 ••• 0^ 




£> = 






1 ••• 








,J = 




• (14b) 








o ... 1) 






= 1,2, 


. . . ,N, are arbitrary complex 






+ ikfT, 


i = 1,2,.. 


, TV, and 





Here aj s ' 1 , s = 1, < 
parameters and i 
ki are complex parameters. We write down the multi- 
soliton solution of the 2-CNLS system^ as below: 
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q s {X,T) = »—, s = l,2 



where 



A I <f> 
-I B T T 

o a o 



(15) 



(16a) 



Here the matrices A and B are defined as 



A 



>j 



ki -\~ 



B 



ki 



i,j = 1,2,..., N. (16b) 



In equation (|16b[) . f represents the transpose conjugate. 
In particular the one-soliton solution (TV = 1 case) reads 
as 



(17) 



q 3 {X,T) 
where j = 1 , 2 and 



An 1 


p-71 




-1 B u 
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. (1) 12 



(fcl + fci*) 2 



(18) 



For the TV = 2 case, we can write the two-soliton solution 
from Eq. (fPoj) as 



Qj(X,T) 



An 


A 12 


1 


e m 


A 2 i 


A22 
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e m 


-1 





Bu B 2 i 








-1 


B\2 B22 


















All 


A 12 1 






A21 


A 22 1 






-1 


Bn B21 









— 1 £?12 £>22 





(19) 



1+ J2 (Pi™e 7n+,f A + fce^+m+n 



Z,m— 1 



(4 



(fc 1 + fc;)(fc 2 + fc;) 

|fcj - fc 2 | 2 (KnK 2 2 ~ Wl2«2l) 
(fc 2 + fc 2 *)(fc 2 +^)|fcl+fc 2 *| 2 ' 

Kim _ \ -y 

~ k + k* 



j) k 2 i - a { 2 j) Ku 1 



'1'/. 



E 



^ + fc* 



(21) 



Similarly, the three-soliton and TV-soliton solutions can 
be written down explicitly. 

Now using the transformation ([3]), the bright-bright 
TV-soliton solution of the two-coupled GP equations can 
be written as 



&j(x,t) =2r V T^exp 



i{6 + fijt) 
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dt 



(22) 



where X and T are given in Eqs. (fl2"b) and (fl2"B). Note 
that the variables X and T are complicated functions of 
the original independent variables x and t as given by 
Eqs. (12c) and (12d). Consequently the width which is 
inversely proportional to R{t), see Eq. (12 c), changes as 
a function t and x, even though it is not so in the case of 
the Manakov system. 



IV. ELASTIC AND INTENSITY REDISTRIBUTION 
INTERACTIONS OF BEC BRIGHT-BRIGHT SOLITONS 

Depending on the various forms of the trap potential, 
gain/loss and interatomic interaction satisfying Eq. (|13p . 
we have deduced different kinds of bright-bright soliton 
solutions using the above expression for the soliton solu- 
tions (|2"2"j) . and the transformations ([5]), (|11|) and (|T^1) . In 
the following, we demonstrate soliton solutions for three 
typical choices of trap potentials. In particular, we have 
focussed on the shape changing and elastic interactions 
of two soliton (TV = 2) and three-soliton (TV = 3) solu- 
tions of the two component BEC. The analysis can be 
systematically extended to arbitrary TV, which we do not 
present here for brevity. 



A. Expulsive Potential 

For time independent expulsive parabolic trap poten- 
tial, il 2 (t) — — Oq, where Hq is a constant, we get R(t) = 



■5 



scch(f2o^ + S) from Eq. (fT3|) . The intensity of the wave 
packet | 2 ) is proportional to R(t)e$ 1 ^ dt as seen from 
Eq. ([22]) and the width of the wave packet is inversely pro- 



portional to J R(t). We have constructed different types 
of soliton solutions by suitably tuning the gain term but 
here we have presented only the constant intensity case 
for 7 = £7otanh(f2o£ + S). Fig. [T] shows one such possi- 
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FIG. 1. (color online) Choice of the atomic scattering 
length a 3 (t) and gain term T(t) as a function of time for 
time-independent expulsive harmonic trap potential V(x) — 
—Q.qX 2 /2, Qo = constant. 



ble gain term T(t) = uj±Qq tanh(f2oi + 5) and the atomic 
scattering length a s (t) = ^asR(t) = \&b sech 2 (fi i + S), 
which can be experimentally realized by varying the ex- 
ternal magnetic field as 



B(t) = B a 



iassech (Q t + 5) 



(23) 



One may note that such a form of scattering length a s (t) 
has been realized in 7 Li and 85 Rb atoms^^. Simi- 
larly the gain term T(t) can be experimentally realized by 
pumping of atoms optically as demonstrated in Refs . 35 ' 37 . 

For N = 1, we get the one-soliton solution of the cou- 
pled GP equation ([XJ) from Eq. (f2"2"j) . As noted above, if 
we choose j(t) = fiotanh(f2of + 5), the intensity of the 
soliton is constant. Figure [21 shows the bright-bright one- 
soliton solution for the two components |1) = |^i| and 
1 2) = for the above gain term where the amplitude 
of the wave packet remains constant. For the N = 2 
case we get the two-soliton solution of ([T]) from Eq. (|2"2"T) . 
In this case, elastic collision occurs only for the specific 



choice of the parameters (a\ ja 



/ (2)/ (2)> 

(a\ ' la\ 



sec 



Eqs. ([2T!]) and (|22p . For all other choices of the parameter 
values, shape changing/matter redistribution interaction 
occurs. Fig. [3] shows the elastic interaction of the bright- 
bright two-soliton solution for k\ = 1 + i, k% = 2 — i, 
a± = = a 2 = c4 = 1- Here the intensity of 



FIG. 2. (color online) One-soliton solution for time- 
independent expulsive harmonic trap potential V(x) = 
— Q.qX 2 /2 and -y(i) = fio tanh(f2o£ + S). The parameters are 
fci = 1 + 0.5i, a[ r> = af ] = 1.0, VL = 890/53 and r = 0.5. 
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FIG. 3. (color online) Elastic interaction of bright-bright two 
solitons for time-independent expulsive harmonic trap po- 
tential V(x) = ~Qox 2 /2 and ~/(t) = fi tanh(£V + S) for 



„(i) 



,.<-) 



,(D 



,.<-) 



A-2 



f2o = 0.06 and ro 



The parameters are k\ = 1 + 
= 0.5. 



the two solitons (Si and S2) in both the components are 
unchanged before and after interaction. The two distinct 
possibilities of the shape changing two-soliton interac- 
tion are shown in Figs.[4ja) and[4|b), see also TableQJa). 
Figs. (Ha) illustrate the shape changing two-soliton inter- 
action for ki = 1 + i, = 2 — i, = a^ 2 ' = = T 
a^ 1 ' = (39 + 80i)/89. Here the intensity of the soliton 
Si gets enhanced while that of soliton S2 is suppressed 
after interaction in the component |1), whereas in the 
component |2) it gets reversed. Fig. HJb) shows another 
possible way of the shape changing two-soliton interac- 
tion for ki = 1 + i, k 2 = 2 - i, a{ 1] = (0.02 + 0.H), 
CV2 = a< i^ = ot% = 1- Here in contrast to the above [cf. 
Figs. 0](a)], the intensity of the soliton Si gets suppressed 
while that of soliton S2 is enhanced after interaction in 
the component |1), whereas in the component |2) it gets 
reversed similar to the well studied case of Manakov sys- 

tem 53 i 55 i 56^ 

For the N = 3 case, we get the three-soliton solutions 
of the coupled GP equation ([T]) from Eq. (|2"2l . Here 



clastic interaction occurs only for a\ 



(i) 



.X 1 ) 



J 1 ) 



.,(2) 



Q 



(2) 



„( 2 ) 



and for all other choice of parameters, 
matter redistribution interaction occurs. Fig. Eta) shows 
the elastic interaction of bright-bright 3-soliton (N = 3 
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FIG. 4. (color online) Shape changing interactions of bright- 
bright two solitons for time-independent expulsive harmonic 
trap potential V(x) = —QqX 2 /2 and j(t) = fiotanh(f2of + S) 
for (a): a' 1 ' = a{ 2) = af } = 1, a£ ] = (39 + 80i)/89 and for 
(b): a[ 1} = (0.02 + 0.1i), = a (2) = a {2) = 1. The other 
parameters are fa = 1 + i, fa = 2 — i, Qo = 0.06 and ro = 0.5. 
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FIG. 5. (color online) Bright-bright three soliton interactions: 
(a) Elastic interaction for = = cup = = 

c4 = c4 = 1 and (b) Shape changing interaction for 
a[ 1} = 0.39 + O.li, a { 2 1} = 0.3 + 0.2i, a { 3 1} = 1.0, a[ 2) = 
(80.0-80«)/89.0, af' = 0.6 + 0.2, a (2) = (89.0 + 80i)/89.0 for 
time-independent expulsive harmonic trap potential V(x) — 
— Q 2 t x 2 /2 and j(t) = f2otanh(f2o£ + <5). The other parameters 
are fa = 1 + 0.5i, fa = 1.5, fa = 2 — 1.5i, fio = 0.06 and 
r = 0.5. 



case) solution for fa = 1 + 0.5«, fa = 1.5, fa = 2 — 1.5i, 
= = «2 = a"^ = 0-3 = a 3 = 1. Here the 
intensities of the three solitons in both the components 



are unchanged before and after interactions. Fig. [5{b) 
shows the shape changing interactions of bright-bright 3- 

soliton solution for af ) = 0.39 + O.li, = 0.3 + 0.2i, 
a 3 1] = 1.0, c4 2) = (80.0 - 80i)/89.0, af* = 0.6 + 0.2, 
c4 2) = (89.0 + 80i)/89.0. Here the intensity of the soli- 
ton Si gets suppressed (S) while that of solitons S2 and 
S3 are enhanced (E) after interaction in the component 
|1), whereas in the component |2) it gets reversed. The 
six distinct possibilities of shape changing interactions 
of three soliton solutions are shown in Table |TJb). The 

TABLE I. Possible combinations of matter redistributions of 
two solitons (N = 2 case) and three solitons (N = 3 case) 
under collision in component |1). E stands for enhancement 
of intensity while S represent suppression. See also Ref.— 
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above three soliton interaction process is equivalent to 
two pairwise interactions. Kanna and Lakshmanan 55 i 56 
have shown that for the corresponding Manakov system 
the first interaction is controlled by the parameters , 
a 2 i a i 2 '; ^1, fa and the second interaction is con- 
trolled by a^, ttg 2 -* and fa. 

Fig. [5] depicts the shape restoring property of soli- 
ton Si during its interaction with the other two soli- 
tons, £2 and S3 for time-independent expulsive harmonic 




t (in ms) t (in ms) 



FIG. 6. (color online) Shape restoring property of Si in three 
solitons interactions for time- independent expulsive harmonic 
trap potential V(x) = — Qqx 2 /2 and 7(f) = ilo tanh(fioi + 5) 
for the choice of the parameters = (39 — 80i)/89, = 
(39 + 80i)/89, 4 1 ' = 0.3 + 0.2i, a? ] = 0.39, a (2) = af = 1, 
fa = 1 + 1.5i, fa = 1.5, fa = 2 - 1.5i, f2 = 0.06 and r = 0.5. 



trap potential V(x) — —£IqX 2 /2 for the choice of the pa- 
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ramcters a 



(i) 



80z)/89, a { 2 ] = (39 



(39 - 

4 2) 



80i)/89, 

(2) 
3 



a K 3 ' = 0.3 + 0.2i, af = 0.39, a\ ' = ap = 1, 
ki = 1 + 1.5i, &2 = 1.5, fc3 = 2 — 1.5i. Here the in- 
tensity of the soliton Si is unchanged after its interac- 
tions with the other two solitons S2 and S3 in both the 
components, while the intensities of solitons S2 and S3 
get changed. The condition for the choice of the pa- 
rameters for the shape restoration of soliton Si is given 
in RcfiSS. Note that the shape restoring property is essen- 
tial for the construction of universal logic gates which are 
necessary for computin g 54 i 56 . The above type of clastic, 
shape changing interactions and shape restoring property 
of solitons during the interactions are similar to the study 
of Kanna and Lakshmana n 55 i 56 in the context of optical 
computing, where the intensities of light pulses are trans- 
formable between two modes. Here, the shape changing 
interactions are interpreted as the transform of the frac- 
tion of atoms between the two components, which can 
be achieved experimentally by suitable tuning of atomic 
scattering length and gain/loss term. Note that exact 
analytical representations for the soliton switching can 
be given in the form of linear fractional tranformations 
leading to logic gates as in the case of optical systems, 
see Ref i 54 i 56 . This type of shape changing soliton inter- 
actions can be used (as disscussed in Refsi2S~— ) in the 
matter wave switching devices, logic gates and quantum 
information processing as in the case of optical comput- 



B. Periodically Modulated Trap Potential 

Next, if we choose the periodically varying atomic scat- 
tering length^ R(t) = 1 + uicos(ujt + 5) and the corre- 
sponding gain term 



7(*) 



ui 2 sm((Jt + S) 
2[l + uJCOs(uJt + S)} , 



u < 1, 



(24) 



with S as a constant, we get R(t) = y/l + u> cos(ujt + S) 
from Eq. ([5]) and the periodically modulated trap fre- 
quency from the integrability condition (| 13[) as 



n 2 (t) 



10w cos( 



3co 2 



1 + cos 2 9 



1 + u> cos 9 



(25) 



where 9 = ojt + S. The intensity of soliton also remains 
constant for this case but the width of the soliton is peri- 
odically varying with time (that is, inversely proportional 

to \/R(t)). Fig. [7] shows the gain term 



r(t) 



lu±lu 2 sin(wt + S) 
2[1 + w cos(ujt + 6)}' 



(26) 



which can can also be experimentally realized by suitable 
optical pumping, and the corresponding choice of atomic 
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t (in ms) 

FIG. 7. (color online) Form of the (a) atomic scattering length 
a s (t) and (b) gain or loss term T(t) as a function of time for 
the periodically modulated harmonic potential with the trap 
frequency (|25[) . 




t (in ms) 



FIG. 8. (color online) One-soliton solution for the periodically 
modulated harmonic potential with trap frequency (|25l) . The 
parameters are fixed as ki = 1 + 0.5i, a[ 
u — 0.4 and ro = 0.5. 



„(2) 



1.0, 



scattering length a s (t) = ^f- [1 + lo cos(ujt + 5)] which 
may be experimentally realized by periodically tuning 
the external magnetic field as 



B{t) = B Q 



A 



— [1 + w cos(u;t + 5)] 



(27) 



Fig. [8] shows the snake-like effect of the one-soliton so- 
lution for the above periodically modulated trap poten- 
tial with R(t) = [1 + uj cos(u;t + d)} and -f(t) as given 
in Eq. (j2~4")l . where the intensity of the soliton remains 
constant in both the components while the width of the 
soliton varies periodically with time. Note that the oscil- 
lation of the soliton goes on increasing with time due to 
the fact that the width of the soliton is inversely propor- 
tional to R(t). It is of importance to note that recently, 
in scalar nonautonomous NLS equation for BEC2£ and 
optical solitons^, a similar kind of snake-like effect has 
been demonstrated. 

Next we analyze different types of two soliton interac- 
tions for the periodically modulated potential for suitable 
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FIG. 9. (color online) Bright-bright two soliton interactions 
for the periodically modulated harmonic potential, (a) Elastic 
collision for = af = = a 2 2 ^ = 1 and (b) Shape 

changing collision for aj 1 ' = aj 2 ' = c4 = lj ~ (39 + 
80i)/89. The other parameters are fixed as fei = I + i, &2 = 
2 — i, uj = 0.4 and ra = 0.5. 



lision of snakc-likc bright-bright three-soliton solutions, 
respectively, for this case. The other parameters are the 
same as in the time independent expulsive harmonic po- 
tential case. The collision effects are similar to the ones 
discussed in the case of expulsive trap potential earlier 
but now the widths of solitons oscillate periodically with 
time because of the periodically varying atomic scatter- 
ing length. 



C. Kink-like Modulated Trap Potential 

Finally, if we choose R(t) = 1 + tanh(cjt + S) and 



7(t) = --[l-tanh(wt + <y)] J 



(28) 



we get R(t) = \J\ + tanh(wt + 8) and the integrability 
condition (fl~3| gives 



n 2 (t) = -i 



i - 



(29) 



(1 +c 2 "') 2 

which is a kink-like modulated trap. For the above case, 




t (in ms) 



FIG. 10. (color online) Bright-bright three soliton interactions 
for the periodically modulated harmonic trap potential. 



Elastic collision for a\ 



(i) 



,.<-'> 



.X 1 ) 



(2) 



J 1 ) 



„( 2 ) 



a 1 = a 2 = a 2 = a 3 



I and (b) Shape changing collision for a{ = 0.39 + 



O.li, a 2 1} = 0.3 + 0.2i, ag' = f.0, a? } = (80.0 - 80i)/89.0, 
a 2 2) = 0.6 + 0.2, af ) = (89.0 + 80i)/89.0. The parameters are 
fci = f + i, ki — f .5 — 0.5i, ki, = 2 — i, u) = 0.4 and ro = 0.5. 



(2) 




fO 

t (in ms) 

FIG. If. (color online) Choice of the (a) atomic scattering 
length a a (t) and (b) gain term F(t) as a function of time 
for the kink-like modulated harmonic potential V(x,t) = 

— LJ 1 — 



(I + e 2 " 4 ) 2 



we sketch the gain 



r(t) = -^[i 



tarm(o;i + 5)] 



(30) 



choice of other parameters. Fig. [^a) shows the elas- 
tic collision of snake like bright-bright two solitons while 
Fig. H^b) shows shape changing collision of snakc-likc 
bright-bright two solitons for u> = 0.4. Figs. [TUTa) and 
(b) depict the elastic collision and shape changing col- 



and the corresponding choice of atomic scattering length 
a s {t) = [1 + tanh(wi + 5)] in Fig. HU The form flU 
of the gain can again be experimentally realized by con- 
tinuously loading the external atoms into the condensate 
by optical pumping as in the case of Refs i 35 i 37 and the 
atomic scattering length can be realized by a kink-like 
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t (in ms) t (in ms) 



FIG. 12. (color online) One-soliton solution for the kink-like 
trap potential. The parameters are ki = 1 + 0.5i, a^' = 
> = 1.0, oj = 0.06 and r = 0.75. 




t (in ms) t (in ms) 

FIG. 13. (color online) Bright-bright two soliton interactions 
for the kink-like modulated harmonic trap potential, (a) 
Intensity unchanged under the collision for at^ = = 
Qfj 1 ' = c4 = 1- (b) Matter redistribution collision for 
c^ 1 ' = a[ 2) = af> = 1, a 2 1] = (39 + 80i)/89. The parameters 
are fci = 1 + i, k 2 — 2 — i, uj = 0.06 and ro = 0.4 




t (in ms) t (in ms) 



FIG. 14. (color online) Bright-bright three soliton interac- 
tions for the kink-like modulated harmonic trap potential, 
(a) Intensity unchanged collision for = = c4 = 
Oj 2 ' — = a ^ = 1- (b) Matter redistribution collision 
for a[ r> = 0.39 + O.li, a 2 1} = 0.3 + 0.2i, q^ 1 ' = 1.0, a[ 2) = 
(80.0 - 80i)/89.0, a { 2 2) = 0.6 + 0.2, a (2) = (89.0 + 80i)/89.0. 
The parameters are k\ = 1 + i,k% = 1.5 — 0.5i, ^3 = 2 — i, 
lj = 0.06 and r = 0.4 



Fig. fT3Ta) shows the bright-bright intensity unchanged 
collision of two solitons while Fig. fT5fb) shows matter 
redistribution collision of two solitons. The intensity 
unchanged and matter redistribution collisions of three- 
soliton solutions for the kink-like modulated potential are 
shown in Figs. H~4T a) and H~4T b). respectively. The other 
parameters are the same as in the time-independent po- 
tential case. Note that, here, the widths of the two and 
three solitons are also decreasing with time. 



V. SUMMARY AND CONCLUSION 



tuning of the external magnetic field as 



B{t) = B 



a° s - ia B [1 + tanh(wt + 6)] 



(31) 



Fig. [T2] shows the one-soliton solution for the kink-like 
modulated trap potential with R(t) = 1 + tanh(w£ + 6) 
and 

7(i) = ~ [1 - tanh(wi + 5)} . 

Here the intensity of solitons in both the components are 
constant while the width is decreasing with time (that is 

inversely proportional to J R(t)). 

Next we analyze the two-soliton (N = 2) and 
three-soliton (N = 3) solutions for the two-component 
BECs with kink-like modulated harmonic trap potential. 



In summary, we have investigated the exact bright- 
bright multi-soliton solutions of the two-component 
BECs with time varying parameters such as trap fre- 
quency, s-wave scattering length and gain/loss term. On 
mapping the two coupled GP equations onto the cou- 
pled NLS equations under certain conditions, we have 
deduced different kinds of bright-bright one-soliton solu- 
tions and interaction of two solitons for time independent 
expulsive harmonic trapping potential, periodically mod- 
ulated trap potential and kink-like modulated potential. 
We have shown the shape changing and elastic interac- 
tions of bright-bright two-soliton and three-soliton solu- 
tions of the two component BECs. The present study 
provides an understanding of the possible mechanism for 
the fraction of atoms transform between the two compo- 
nents. Especially the shape changing collisions of matter 
wave solitons can used for matter wave switching devices, 
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logic gates and quantum information processing. These 
elastic and shape changing interactions can be realized in 
experiments by suitable control of time dependent trap 
parameters, atomic interaction and interaction with ther- 
mal cloud. 
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